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Abstract 

Proceeding from the simplest field theoretical model of 2D metal, the normal phase 
Green functions of the weakly interacting fermions and the order parameter fluctua- 
tions (responsible for the attraction between fermions) are obtained. It is shown that 
taking into consideration the fluctuations mentioned leads to a considerable reduction 
of the fermion wave function renormalization constant (quasiparticle weight) as well 
as to a linear dependence of the quasiparticle damping on the temperature. A general 
dependence of 2D Fermi liquid properties on the fermion density is discussed. The rel- 
evance of the proposed model to the marginal behavior of the Fermi liquid of high-T c 
superconductors, in particular, to their linear temperature dependence of the resistivity 
is indicated. 
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1 Introduction 



Theoretical study of the normal properties of 2D (as well as quasi-2D) metals was strongly 
intensified by the discovery of high-T c superconductors (HTSCs). These compounds display 
a number of nonconventional features that make it difficult to use the Fermi liquid (FL) 
model to describe their main properties. The most striking features are: (i) an essentially 
reduced quasiparticle weight (QW), Z, of the excitations in the vicinity of the Fermi surface; 
(ii) a rather large damping of the corresponding quasiparticles; (iii) a linear dependence of 
the resistivity on the temperature above T c , and many others (see, for example |T|). Now, as 
far as we know, there is no satisfactory understanding as to what makes these compounds 
behave so strangely[] and why they cannot be described by the theory of Landau FL. Thus, 
searching for theoretical models, being able to explain at least some of those odd properties 
of HTSCs, remains a very important problem. In this paper we shall make an attempt to 
show that one of the possible answers is that the observable properties of HTSCs can be the 
result of the extreme susceptibility of 2D attractive interaction FL with respect to a deep 
intrinsic structural rearrangement. 

We recall that the most characteristic property of the ideal Fermi gas is the existence 
of the Fermi energy, ep, which, in accordance with the Pauli principle, separates occupied 
and empty states. The distribution of particles over the energy (and momenta) at T = 
can be represented by the step-function. Adding to the system a weak interaction between 
the particles corresponds to the transition from the ideal Fermi gas to the FL. When the 
interaction is weak and there is no danger of any kind of instability, the distribution function 
changes only a little bit: the QW slightly decreases (Z ~ 1) and a weak damping, 7 ~ 
(e — e^) 2 , appears |]2], [J| . In other words, despite the interaction between fermions, the Fermi 
surface, determined by the "jump" in the distribution function, does exist in the system 
(this is the direct manifestation of the famous Migdal theorem). 

The situation changes drastically at lower dimensions and/or in the case of a very strong 
interaction between fermions. As for the low dimensions, the most convincing example is the 
ID conductor, where even the weakest electron-electron interaction causes the disappearance 
of the jump in the distribution function j|. Such a FL (called so Luttinger), strictly speaking 
cannot be referred to as FL since all its excitations belong to those of the bose-type [fflBI, BJ . 

The example, demonstrating the influence of an interaction, is the strongly correlated 
Hubbard model. On-site repulsion, U, in this model acts on all particles in the system in 
contrast to the case of a generally accepted FL where the interaction is actual only in the 
small vicinity of the Fermi surface. As to the result, such a repulsion leads to a considerable 
reduction of the QW. When the repulsion increases — the QW decreases 0, [§| and in the 
limiting case U — > 00, the Fermi surface disappears completely [Z = 0) at any fixed density 
of fermions || |10j. 

As for the HTSCs, according to much experimental data, neither Landau no Luttinger 
FL can explicate them. What is interpreted as QW in the case of HTSCs is relatively small 



(Z ~ 0.4 — 0.6 [II], 0]), so the corresponding FL could not be considered as Landau FL. 



1 Sometimes they are called the "strange metals" 
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On the other hand, it is not Luttinger FL since the QW is not small enough. In addition, 
quasiparticles themselves are less stable than in the case of ordinary Landau FL. Actually, 
the damping of them in a close vicinity of the Fermi surface at T = looks like 7 ~ |e — €f\°, 
with a ~ 1 |12], [13| ; this differs from the usual FL power law with a = 2. Such a FL with 
these exceptional characteristics and a qualitatively different "intermediate" behavior is now 
called "marginal FL" |§ [Lg ||, 

If one assumes that the marginal FL is relevant to HTSCs, then a sufficient number of 
observable physical properties of HTSCs can be quite naturally explained. Up to the present, 
however, the latter has been achieved only in the framework of some phenomenological theory 



14| , containing a lot of additional assumptions, the physical meaning of which has yet to be 
clarified. 

In what follows we consider the simplest 2D field theoretical model in the normal phase 
with a weak inter-fermion attraction (its nature is not specified here). As will be seen, 
the Green function of the order parameter fluctuations (OPF), responsible for the attraction 
between fermions in our model, turns out to be consistent with the very general requirements 
of the phenomenological theory |TJ|, being capable of explaining many HTSC properties. 
Here we should note, however, that in the cited paper the interaction between fermions 
put down due to the exchange by charge or spin fluctuations. Though it is crucial for the 
interpretation of some experimental data, as well as for the calculation of the Green function 
itself, the latter does not play any important role. Therefore, below we shall restrict ourselves 
to studying only those properties of HTSCs which can be extracted from the Green function 



of interacting fermions. In order to be able to consider the other properties, as in Ref. [[14 
we would need to calculate corresponding charge and spin correlation functions. But that, 
in fact, is beyond the scope of this paper. 

Below it will be shown that taking into account the interaction mediated by the OPF leads 
to a considerable reduction of the QW even at a very weak attraction between fermions. This, 
we believe, is the result of dimensionality lowering, what is accompanied by the enhancement 
of the role of any attractive interaction. Investigating the dependence of the QW on the 
density of fermions, we shall come to a quite natural conclusion that decreasing the density 
leads to decreasing the QW and increasing the damping. Note that just these two features: 
(i) a quasi-2D character and (ii) a relatively low density of carriers, are as a rule emphasized 
when properties of HTSCs are discussed |jlj]. 

Lastly, maybe the most interesting result of our paper is a linear temperature dependence 
of quasiparticle damping. It is just this feature that was in fact a priori accepted in the 
phenomenological description \T%\ for explaining the surprize linear temperature dependence 
of resistivity observed in high-T c compounds. 

Throughout the paper we use units in which h = 1 and fee = 1. 
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2 Model and general discussion 



To solve the problem analytically as far as possible, we choose the simplest field theoretical 
model of fermions with attraction. The corresponding Hamiltonian is ||: 



H = J d 2 rH(r), 



(1) 
(2) 



Here we use the standard notations: m — for the effective fermion mass, \x — for the chemical 
potential, o —\ or \ denotes a spin variable. The interacting constant, g, is positive what 
corresponds to the attraction between fermions. 

Introducing Nambu spinors for fermion field ^ = (ipt,ipi), one can rewrite Eq.(|2|) in a 
more convenient form: 



H(r) 



(3) 



where r z , r + = (t x + ir y )/2, r_ = (r x — ir y )/2 are Pauli matrices. 

Below we shall use the functional integral approach along with the Matsubara thermal 
technique. Thus, the partition function is expressed through the Hamiltonian as: 



Z(T) = J [d&dV] exp - Jdr J d 2 r (V<9 T * + H(r] 



(4) 



where T is the temperature and [d^'d^] denotes the measure of the functional integration 
over the Grassmann variables \& and W, satisfying the antiperiodic boundary conditions: 
\tr( r; r ) = -Hf(r + (3; r) and ^(r; r) = -^ f (r + /3; r). 

To be able to construct arbitrary Green functions one should add to (f|) classical sources 
associated with the Grassmann fields in the Hamiltonian (Q). However, we shall not write 
them down explicitly, even though it would facilitate the understanding of some key formulae 
below. 

Using now an auxiliary Hubbard-Stratonovich complex scalar field in the usual way, one 
can represent the Eq.(H) in an equivalent form: 



Z(T) = J [d&d$d$d$*] exp - J dr J d 



l$l 



+ 



(5) 



The main virtue of this representation is a nonperturbative introduction of the composite 
fields: $ = —g^^T+ty and $* = — gWr+ty and possibility to develop the consistent ap- 
proach. Specifically, the expression (|5|) turns out to be rather convenient for studying such 
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a nonperturbative phenomenon as, for example, superconductivity. In this case a complex 
Hubbard-Stratonovich field naturally describes the order parameter arrising due to forma- 
tion of pairs (to some extent, the Cooper ones), or composite bosons. The average value of 
|$| is proportional to the density of these bosons, on one hand, and determines the gap in 
the one-particle fermi-spectrum, on the other. In the normal phase the homogeneous order 
parameter disappears, i.e. |$| equals zero, but the OPF inevitably take place even there. 
Intuition, however, suggests that these fluctuations should play the crucial role in a small 
neighbourhood of the corresponding phase transition. As it will be shown, this is in fact 
the case in our model. Moreover, OPF not only exist but also have an influence on different 
characteristics of normal phase on the whole. 

Before going further, it may help to recall the most essential property of 2D systems. 
The question is the famous Mermin-Wagner-Hohenberg (MWH) theorem which forbids the 
appearance of a complex (in our case — charged) order parameter $ at any finite temperature 
|T7| , pi| . If we studied the problem in the mean field approximation we would obtain a finite 



critical temperature T c , for instance, of the superconducting phase transition. Taking into 
account the next to leading order approximation immediately reduces T c to zero. Such a 
manifestation of the MWH theorem may make the study of pure 2D superconductors at 
T^O meaningless. However, as was shown recently, the situation here is much more subtle. 
Actually, 2D metals undergo two phase transitions [jnj. The first one is connected with 



violating only some discrete symmetry and the appearance of a neutral order parameter 
(modulus of the Hubbard-Stratonovich field, |$| ). It is this transition that was usually 
interpreted in the mean field approximation as the appearance of superconductivity. Despite 
this modulus being different from zero, the field $ on the average is zero, due to the temporal 
and spatial fluctuations of its phase factor. We would call such a state of the system by the 
term "molecular gas", emphasizing by this notion that there is no phase coherence in the 
"Cooper pair" dynamics^. The second phase transition is connected with the appearance of 



the so called Berezinsky-Kosterlitz-Thouless (BKT) phase pffil . Developing of such a 
phase is just accompanied also by the superconductivity |T9|j . 

In what follows we shall study only the normal phase of 2D metal, i. e. the phase where 
even the neutral homogeneous order parameter is absent. We believe that precisely this phase 
of pure 2D model represents the true normal phase of quasi-2D HTSCs more adequately. 

Below we shall use a generally accepted assumption that the approximation taking into 
account only the quadratic (Gauss) fluctuations of <3> (around zero in our normal phase) 
describes the system quite well. Thus, integrating out the fermion fields in ([|), one obtains: 

P P 



Z(T) = exp{-S ) J[d$d$*]exp - J dr x J dr 2 



o o 



x J d 2 r J dVr^rOr 1 {t x -r 2 ;ri - r 2 )<£(r 2 ; r 2 ) 



(6) 



2 It must be noted that in Ref. |2(| there was introduced the temperature T p which characterized the 
energy scale of composite boson formation. 
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where 



is the one-loop effective action, and 



So = TrLn[G 



(7) 



r _1 (r;r) = -5(r)5(r)+tr 
9 



G 



0\T 



r)r_G (-r; -r)r + 



is the inverse Green function (propagator) for OPF. Both quantities (|7|) and 
through the free fermion Green function which obeys the equation: 



- d T + t z 



2m 



G (r;r) = 5(r)5(r) 



with the following boundary condition: 

G (r + /3; r) = -G (r; r). 
The integration in can be formally performed. The result is 

Z{T) = exp(-S - TrLn[Y- 1 }). 



are expressed 
(9) 

(10) 
(11) 



If we could have obtained the explicit expression for the partition function we could, in 
principle, obtain all thermodynamical functions but in reality it is usually impossible to do 
this. Therefore, we shall study more simple objects, namely, the Green functions which can 
also give the information about the system. 

To avoid possible further misunderstanding, we shall write down the formulae which are 
used throughout the paper for Fourier transformations, connecting coordinate and momen- 
tum representations: 



< 

F(iu n ,\i) = J " dr J d 2 rF(r, r) exp(zcj n r — zkr), 
o 

. . m ~ r d 2 k 
F(r,r) = T > / — — F{iu n , k) exp(-zuj n T + ikr) 



(12) 

(13) 



where u n = %T{2n+ 1) are fermion (odd) Matsubara frequencies. In the case of bosons they 
should be replaced by even ones: Q n = 2nTn. 



3 Green function of OPF 

The Green function for free fermions in the momentum representation (|12|) takes, as follows 
from (|9|), the very simple form: 

G (iu n , k) = l j- 2 r- (14) 
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Substitution of this expression into (BI) leads to the explicit formula for the inverse Green 
function for OPF in the momentum representation: 

T,-i/,n 1 1 r «,tanh(^+/2)+tanh(^-/2) 



zfi n ; K) = - - — / d 2 * ~ 7 ' (15) 

where £±. = ^(kiK/2) 2 — //; k and K are relative and total momenta of a pair, respectively. 
The divergence at large momenta, k, in (|T5|) is regularized by the ultraviolet cutoff, W, which 
may be considered as the free fermion band width. Assuming that the cutoff is much bigger 
than all other characteristic scales of the problem, we can get rid of it after coupling constant 
renormalization and next turning W to infinity (see, for example, p3| ). So, expressing the 
coupling constant in the form 

1 m , 2W 

- = 7 -ln TT , 16 
g An \e b \ 

where e& is the two-particle bound state energy, and substituting it into Eq.flTSp, one can 
immediately take the limit W — > 00. The detailed calculation can be found in Appendix 
A (formula fliT])). As is easily seen, the expression obtained is rather complicated for use 
or even for qualitative analysis. However, there is one particular aspect of it which should 
be emphasized here — the asymptotics of the imaginary part of T~ 1 (fi;K). The explicit 
expression for it is: 

Imr -1 (fi;K) = -— ^(fi + 2/i-K 2 /4m)tanh(/3fi/4) x 

x[l - f(n + 2i2-K 2 /Am,K)}, (17) 

where the function f(x, y) (see Appendix A) varies from zero to one and vanishes as x, y — > 
or x, y — > 00. It is easy to be convinced that the asymptotics of ( P7|) at Q — > and Q — > 00 



coincide exactly with those postulated in Ref. [14]. Recall that the phenomenological theory 



which was proposed in that paper explains in principle a lot of experimental data on HTSCs. 
Therefore, one can hope that our microscopic approach will also be able to describe correctly 
some observable properties of HTSCs. 

In order to calculate the Green function of interacting fermions (see the next section) 
one still needs to find some rather simple approximation for OPF Green function. The 
most commonly used one is the so called the derivative expansion, which is valid at small 
frequencies and momenta. To obtain a consistent expression we need first to construct the 
retarded real-time Green function (by means of analytical continuation ifl n — > Q + iO) and 
then to expand it at K — > and Q — > 0. Using the explicit formula ([41]), we arrive at the 
following representation: 



K 2 

- 1 (n;K) = a + b dl (1<< 

K ' Am y 



where 



rn 
Arc 



, 7T r dx x(3a 
In — — : / — tanh 



/3|e 6 |7 J x 2 







(19) 
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87T/i 



in 



tanh^ + ^ + ^ 
2 2 4 



dx 



tanh x 



-/9m/2 



^tanh^ + A(2 M ) 



.1T 



(20) 
(21) 



and In 7 ~ 0.577 is the Euler constant. Note that here and henceforth the same symbols are 
used for both the real-time and thermal Green functions. 

The remarkable feature of the expansion fll8[) is the existence of a term that is linear in 
frequency (the next, quadratic, term is omitted in flT8|)). We stress this point here because 
sometimes the papers devoted to nonrelativistic systems appear (for example, p4| ) where 
even the possibility of such a term in derivative expansion (or, what is the same, in the 
composite boson effective action) is not assumed. Perhaps this is the result of the uncritical 
use of some relativistic methods in nonrelativistic systems. In any case, there is no any 
reason for omitting this term in the model under consideration. Note that the importance 
of such a term has been recently emphasized also in Ref. P?5[] , where a similar model 
for 3D metal was investigated. Finally, the appearance of a term that is linear in frequency 
seems to be quite natural in Galilean-invariant models, and its absence should be regarded 
as being suspicious, if there is no any rather convincing reason for its absence. 

We note here that the condition, indicating the appearance of instability in the system, 
reads as a(/x, T) = 0, i.e. (see (|T9"D) in our caseQ: 



In 



7T 



/3|e 6 |7 



dx xBil 

— tanh 

x 2 



0. 



(22) 



This equality defines a curve (below it will be called "the stability line") which separates 
two different regions on the /i-T plane (see Fig.l). The region "N" in Fig. 1 corresponds 
to the normal phase of the system where there is no order parameter. Inside the second 
region, denoted by "OP", some nonzero order parameter appears. It has to be noted that 
the developing of the BKT phase also takes place somewhere inside the OP area. 

As for the explicit dependence \i = /x(T) on the stability line determined by the equation 
([22]), we can only give asymptotics at both low and high temperatures: 



2 

7T 2 T 2 



at 
at 



T> leJ 



(23) 
(24) 



To understand the physical meaning of the stability line one needs to depict the curve 
which describes the physical dependence of the chemical potential (so far it was treated as 

3 In Ref. |^3| just such an equation was considered as one of the system of self-consistent equations for 
determining the values of T c and fi c = fJ*(T c ) in mean-field approximation at constant density of particles. 
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an independent external parameter) on the temperature at fixed density of carriers in our 
system. The equation that determines this second curve is 

n A^ T ) = -nryl / d 2 p tr[r z G(iuj n ,p)], (25) 

V / n=— oo 

where G(iuj n , p) is the interacting fermion Green function expressed in its turn through the 
temperature and the chemical potential. Note that, after taking into account the expression 
(PD and the definition for the bare density: Uf — T/Vd In Z(T) / dfi, Eq. ( pop can also be rep- 
resented as the sum of two terms (see also |25|, ^Tj). The first term, rip, is expressed through 



the free fermion Green function and the second one, ub, — through the Green function of 
OPF (or composite bosons). Thus, one may interpret such a situation as if a dynamical 
balance rif({i,T) = np(/x, T) + 2ne(/z, T) in the distribution of the initial noninteracting 
particles over fermion and boson degrees of freedom is established in the system. 

As the first approximation, we could substitute the free fermion Green function (^) into 
(f25p what is equivalent to omitting the bosonic term in the expression for the distinctive 
"conservation law" rif = rip + 2ns- As a result, one can see that the curve which describes 
the dependence of the chemical potential on the temperature looks somewhat similar to the 
dashed line in Fig. 1. Increasing or decreasing the density of carriers in the system would 
lead to moving this curve up or down, respectively. 

Turning back to the physical meaning of the stability line, one can be convinced that, 
changing the temperature in the system corresponds to removal to the right along the dashed 
line in Fig. 1. Thus, we can naturally conclude that our system is in its normal phase at 
temperatures that are high enough. On the contrary, decreasing the temperature corresponds 
to displacement to the left along the dashed line. When the system proves to be in the point 
where this line intersects the stability line a neutral order parameter appears. Subsequent 
removal, in OP region, will result in the transition into BKT phase (the second stability line, 
which must be calculated from another equation lfil|, is not depicted in Fig. 1). 



4 Fermion Green function 



The general form for the inverse thermal Green function of interacting fermions reads (com- 
pare with fll4j)): 



G l (iu n ,p) = iu n - t z 



2m 



(26) 



where the self-energy operator for interacting fermions, H(iu n , p), can be already expressed 
through the free Green function of fermions and OPF Green function by means of the 
following formula: 

d 2 k r 

r + G (iu} n >, k)r_r(zO n _ n /, p - k) + 



p) 



7~z \ ^ 
J J- 



n'=— -oo 



(27T) 2 l 

+T-G (iuj n r, k)r + r(zf2 r 



k- P ) 



(27) 
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To trace the complete derivation of this expression one needs to use the standard Feynman 



technique for the system given by (g) |28| . Now the substitution of the OPF Green function 
in its derivative approximation flT8|) into the last formula leads to the integral representation 
for the self-energy operator (f27|): 

where the summation over the Matsubara frequencies was performed. Assuming that the 
result of integration in the last expression does not strongly depend on the particular form of 
the smooth function in the numerator, we can obtain the following approximate self-energy 
operator: 

E(<un, P) = 7 (2n)2 (p _ k)a '(J™' S. . • (29) 



To arrive at this expression, we have carried out some formal limit of zero temperature in (|28|). 
In other words, only the explicit dependence of this expression on the temperature through 
the arguments of the hyperbolic functions was subjected to this limit: coth /?(...) — ► 1, 
tanh /?(...) — > sign(...). Of course, such a treatment introduces some uncertainty in the 
quantitative analysis. But the advantage of making further step by means of analytical 
manipulations is so attractive that we would like to expect the reliability of such an approx- 
imation. After simplification mentioned, the integration over momenta in Eq . (|29|) can be 
performed explicitly and the result is written down in Appendix B (formula (|44|)). 

The real-time retarded fermion Green function in a close vicinity of its pole (as follows 
from (|26|)) has the standard form: 

G(w,p) = t| (30) 

where Z is the introduced above QW, 7 P is the quasiparticle damping and E p is some function 
determining the low-energy one-particle spectrum. All these functions in ( p0| ) are expressed 
through the self-energy operator (p9[ ): 



1 - r z Re , (31) 




E p = Z[l + 2mRe ^iy j ) [ ) , (32) 

7p = -ZIm ^S(0;p ) + "^"V ~Po) j , (33) 

where po is renormalized "Fermi" momentum, which is the solution to the equation E p = 0, 
or 

A = tl - Re x(0; Po ). (34) 
2m 
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Note that this equation could be solved, in principle, at any chosen values of the chemical 
potential and the temperature. However, taking into account physical reasons, we restrict 
ourselves to solving the equation (|34T) for two sets of parameters {T, /i} only. The first one 
corresponds to those values of the temperature and the chemical potential which lie on the 
stability line. From the physical point of view, this choice implies that the system is kept 
close to the critical region, but the density of carriers in the system increases monotonically 
at moving up along the stability line, or decreases in the opposite case. As a result, one 
will obtain the qualitative dependence of the near-critical properties of 2D metal with inter- 
fermion attraction on the density of carriers. The second set relates to the curve (the dashed 
line in Fig. 1), where the density of carriers is constant. This gives us the information about 
the dependence of 2D metal properties on the temperature. 



5 Numerical results and discussion 

We start from the analysis of near-critical properties at different carrier densities. As was 
noted above (see Section 3), the stability line is determined by the condition a(/i,T) = 0. It 
is easy to check that being taken into account this condition simplifies the analysis essentially. 
Indeed, the only quantities which one needs to know for the calculation of Z, E p and 7 P0 
(defined in (|30|)) are the self-energy operator and its first derivative with respect to frequency 
at uo = (see (|31|)-(j33]) ). Both these quantities are given in Appendix B. Performing the 
analysis, we come to the conclusion that in the case of parameters belonging to the stability 
line the relevant self-energy operator does not depend on the momentum and, as a result, 
the equation (|34]) transforms into expression for the momentum p . Substituting then this 
momentum into (|3T|), fl3"2|) and ( |33| ) we would obtain complete information about the Green 



function for the interacting fermions. We, however, shall not write down corresponding 
cumbersome complex expressions, but instead present the graphs for £ = (po 2 /2m — fi) /fi, 
Z and 7 Po as functions of (3[i (Fig. 2-Fig. 4, respectively). 

Note that only the region of positive values of the chemical potential is considered. 
Our investigation fails in the region near the point \i = 0, since the derivative expansion 
(|T8|) is not valid there. Though there are no obstacles for studying the region of negative 
values of the chemical potential separately, we shall not do that because we are interested in 
obtaining some qualitative information about the properties of the system proceeding just 
from the general monotonic dependence of the chemical potential on the particle density on 
the stability line. The singularity in Eq . (|18|) at \i — (because of c — > oo, see fl2~T[) ) forbids 
us to go into the region of negative chemical potential continuously. 

Since the parameter /3fi monotonically changes along the stability line, the graphs in 
Fig. 2-Fig. 4 can be thought of as a qualitative dependence of corresponding quantities on 
the density of carriers in the system. Thus, we conclude that when the carrier density is 
high the system naturally tends to become conventional FL with a rather small damping of 
the particles and QW practically equals one. In the opposite case, damping increases and 
QW reduces, i. e. marginal features of FL really appear. 

Let us now turn to studying the properties of the system at fixed density of carriers. This 
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case turns out to be much more difficult than the previous one. In particular, the equation 
([34]) cannot be solved so easily as before. In fact, one should consider the system of two self- 
consistent equations fl25|) and (|34]) at each chosen temperature. To simplify the system, one 
can use a quite natural approximation in which the additional equation (^) is taken in the 
mean field form, i.e. instead of interacting fermion Green function we shall use the free one. 



But even in this case the equation (|34j) cannot be solved explicitly. Solving it numerically 
and substituting into (|3~1~D and (|3~3|), we obtained the dependence of £ = (p 2 /2m — /i) //i, 
Z and 7 po on the temperature. The corresponding graphs are represented in Fig. 5-Fig. 7 
for several different densities: (1) ep/\eb\ = 500; (2) e.p/|e&| = 750; (3) / 1 e& | = 1000, 
where €f = irn/m by definition. All curves start at the corresponding critical temperatures. 
For example, it is seen that in the wide range of T's the QW considerably differs from 1, 
essentially depending on T. 

But the most amazing, to our mind, result of the numerical analysis is a linear temper- 
ature dependence of the quasiparticle damping within a rather wide range of temperatures 
(from T c up to about 4T C and even higher). As was noted above, it is just such a linear de- 
pendence that was supposed in Ref. |TJI] for explaining the unusual temperature dependence 
of resistivity of HTSCs^]. Our calculations show also that the decreasing of the ratio €p/\€b\ 
rather quickly breaks the linear dependence of damping 7 po on T what can be probably 
explained by the sharp decreasing /i up to so small values where the derivative expansion 
fails. 

Of course, numerical results presented here cannot, by themselves, give us a completely 
convincing answer concerning the physical reasons for an explanation of the peculiarities of 
the 2D metal model. Our hypothesis is quite simple and seems to be rather natural: the 
most important features which are responsible for the marginal FL behavior are lowered 
dimensionality and, of course, attractive character of inter-particle interaction. Indeed, the 
two-dimensionality has resulted in the asymptotics of the imaginary part of the inverse 
Green function (|l7l): 

Jmr-VtyK) ~ at ( 35 ) 

16 

ImT-^K) ~ - — , at Q -> oo, (36) 

which were required in Ref. [0 (pay attantion that the Green functions for spin and charge 
fluctuations were used in the paper cited). As to the attraction, it may play the important 
role even in 3D models but in 2D case its effect proves to be decisive. While in the former 
case OPF are very important mainly in the superconducting phase, in the latter one they 
turn out to be of great importance in the normal phase too. 



4 This result does not seem to correspond to the marginal Fermi liquid definition (see Introduction). It 
must be kept in mind, however, that the definition mentioned refers to the case T = which should be 
excluded from our consideration by physical reasons. As to T ^ the quasiparticle damping proves to be 
nonzero (see (|33|)) even at Green function pole p = po- 
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6 Conclusion 



In this paper we studied a model of a 2D metal with a local weak attraction between fermions 
in the normal phase. Within the framework of rather natural approximations, we have 
shown that such a system displays properties which differ essentially (mainly in the vicinity 
of the stability line) from the conventional properties of Landau FL. In particular, the QW 
noticeably differs from 1 and the damping of quasiparticles is rather large. The manifestation 
of this effect becomes more essential as the density of particles decreases. 

We also investigated the temperature dependence of the QW and the quasiparticle damp- 
ing. The most important result, we believe, is the behavior of obtained curves for the quasi- 
particle damping, exhibiting features which can explain the linear temperature dependence 
of resistivity which is observed at all temperatures of HTCS normal phase existance. 

Since through the paper we restricted our consideration of the 2D model with non- 
retarded inter-particle interaction only which, as is well known, is valid in the case of rather 
small Fermi energy, it seems to be very interesting and important to investigate the influence 
of an indirect interaction between fermions (mediated, for example, by the phonons, magnons 
or some other bosons and resulted in a retarded attraction) and its possible connection with 
the OPF in the normal phase too. One more rather urgent for HTCSs problem is to study 
the possible effects when the OPF have non-trivial (for example, c?-wave) symmetry. At 
last, we would like also to note that since our approach involves some intuitive assumptions 
(even though they are commonly used in solid state physics), more consistent consideration 
would be useful. 

One of us (V.P.G.) would like to acknowledge Prof. D.Atkinson for critical reading the 
manuscript. The research was supported in part by the Grant No.K5O100 of the Inter- 
national Science Foundation and by the International Soros Science Education Program 
(ISSEP) through the Grant No.PSU052143. 

Appendix A 

Here we derive the renormalized OPF Green function. Substituting the coupling constant 
(|n5|) into flTSP, one obtains 



r-^o^K) 



m r, 2W 



2W 

J dx 
o 



tanh ( ) 9( g + g-2/x)/4) 




+ 



47T . 




(37) 



where the function f(x, K) is 




cosh (p(x + £ - 2//)/2) + cosh ( 



cosh(f^cos0)-l 



COS0 



(38) 



71 J 
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(see also (|i~T|)). As was pointed above and can be easily checked directly this function varies 
from zero to one and vanishes when x or y — > as well as when x or y — > oo. 
Using identities: 



1 



x + ^ - 2/i - iVt n 



4 



+ 



+ £-2*0 (^ + £-2/i-zO r 



.r 



(39) 



and 



— tanh x = In 

X 7T 



/ — (tanhx — 1), at a > 0, 
7 a; 



(as earlier In 7 is the Euler constant), one can take the limit W — > 00 in 
the final expression: 



(40) 

and come to 



777 

r-^jK) = — 

47T 



7T f dx ( K 2 

In — — + / — tanh x(3(— - 2/i)/4 
p|efe|7 J x \ Am 



-in. 



+ 



rfxtanh (p{x + g - 2/i)/4 

«( (* + £-2/x) (^ + £-2/^-^ 

cfetanh + g - 2/i)/4 



+ 



ac + ^ - 2/i - zfi„. 



(41) 



Appendix B 

To calculate the integral (p9|), one needs to use the following table integrals (see, for example, 

at \b\ > \a\, (42) 



7T 



and 



b 4- a cos <p ^Jb 2 — a? 



dx 1 



R(x) 



\n(2JcR(x) +2cx + b , 



(43) 



where i?(x) = a + 6x + cx 2 and c > 0. 

Thus, straightforward manipulations lead to the formal result for the self-energy operator 
expressed through the elementary functions: 



b — c 
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where the dimensionless functions 



Atc ~ Air \x Att /J> 

a((3\eb\, /3[i) = — a, b(j3fi) = b, c{(3^) = c, (45) 

m mm 

S = -S, £=l(lf--r) t Q=^u, (46) 
\i fj, \2m J l\i 

more convenient for numerical calculations, were introduced. Note that one of the parameters 
(c) in (fHf ) is complex, so there is some uncertainty in the expression. To avoid it one needs 
to clarify which branch of the complex function is chosen. However, in the main text of the 
paper only the self-energy operator and its derivative with respect to frequency at uo = 
and large value /3// are used. Both required functions at those conditions are free of the 
mentioned uncertainty, so we shall not concern about the difficulty of choosing the needed 
branch. 

When parameters of the system are chosen so that they obey the equation a(/3|e&|, = 
0, which identically coincides with ( p2|) and determines the stability line, the self-energy 
operator and its derivative at ui = become very simple: 

S(o ; n = ^-^hi- at |<0; (47) 
b-c c 



dr z uj 



Ac 



4b 



uj=0 



fKi + c) (b-c)(£b + c) 



at 



£<0. 



(4* 



The expressions for £ > can also be obtained from (fPj). We do not write down them 
since they were not used in our analysis in the main part of the paper because the equation 
£ = — S (see Rl) only has a solution with negative sign. 
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